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Abstract 

In bilayer systems electron-hole (e — h) pairs with spatially separated 
components (i.e., with electrons in one layer and holes in the other) can 
be condensed to a superfluid state when the temperature is lowered. This 
article deals with the influence of randomly distributed inhomogeneities 
on the superfluid properties of such bilayer systems in a strong perpendic- 
ular magnetic field. Ionized impurities and roughenings of the conducting 
layers are shown to decrease the superfluid current density of the e — h 
pairs. When the interlayer distance is smaller than or close to the mag- 
netic length, the fluctuations of the interlayer distance considerably reduce 
the superfluid transition temperature. 
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1 Introduction 

The bilayer systems with pairing of electrons of one layer with holes of the 
other layer have attracted considerable interest [TJ [H [3] because such electron- 
hole (e — h) pairs can be condensed to a superfluid state, accompanied by a 
counternow superconductivity (as predicted in Refs. 4,5). This special type of 
superconductivity manifests itself as equal and antiparallel nondissipative cur- 
rents flowing in the conducting layers when the e — h pairs move in one direction. 
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Experiments on antiparallel currents were carried out with GaAs/AlGaAs het- 
erostructures formed by two closely spaced quantum wells, each containing a 
quasi- two-dimensional electron gas [5J [7J |S]. It has been found that in a per- 
pendicular magnetic field, when the total Landau level filling factor ut is equal 
to unity (vt = v\ + ^2 = 1), a- lowering of the temperature leads to a signifi- 
cant (more than an order of magnitude) decrease of the longitudinal resistance 
p xx and the Hall resistance p xy in each layer. The obtained results can be in- 
terpreted in terms of superfluidity of the e — h pairs with spatially separated 
components (interwell excitons) [7J Indeed, the condition Vt = 1 is equiva- 
lent to the equality v\ = \ — vi which means that the number of occupied states 
in the first layer exactly coincides with the number of vacant states (holes) in 
the second layer. If the temperature T = 0, all electrons from the first layer 
and all holes from the second layer form e — h pairs and there are no unpaired 
carriers. Due to superfluidity of the e — h pairs the longitudinal resistance p xx 
in each layer is equal to zero. The Hall resistance p xy must also become zero 
because e — h pairs (being electrically neutral) does not contribute to the Hall 
effect. The residual values of p xx and p xy observed in the experiment (at T^O) 
are caused by motion of free vortices and unpaired carriers respectively. 

The properties of the superfluid state of e — h pairs in the bilayer system 
depend critically on the interlayer distance d. Experiments show that the de- 
crease of p xx and p xy described above is not observed if d exceeds I.&Ih - 1.91h, 
where Ih = \Jhc/eH is the magnetic length [6l [8]. According to theoretical 
concepts [10] , with increasing d the coherence of e — ft pairs breaks down if the 
coherence length becomes smaller than the average distance between the pairs. 
This transition is of a quantum nature, since at large d the interlayer coherence 
is destroyed by quantum fluctuations 

The main attention in this article is paid to the properties of systems with 
relatively small interlayer distances d < lu- This range of interlayer distances 
is quite important because in a homogeneous system the critical temperature 
of the superfluid transition (T c ) reaches its maximum value at d rs OAIh (see 
later). The maximum is rather sharp, since with reducing d the interaction 
between the e — h pairs is weakened, and the value of T c decreases rapidly due 
to thermal fluctuations. It turns out that in an inhomogeneous bilayer system 
with small d the weakness of the interaction between the e — h pairs leads to a 
high sensitivity of T c to the structure inhomogeneities. As a consequence, the 
critical temperature is strongly lowered by inhomogeneities for d < Ih- (See 
Ref. [12] and the references cited therein for the influence of inhomogeneities 
on the bilayer system at d 3> Ih-) 

2 Suppression of superfluidity of e — h pairs by 
the electric field of ionized donors 

In the bilayer systems that are based on semiconductor heterostructures, elec- 
trons fill of the quantum wells due to ionization of donor atoms. We will show 
that the electric field of the donor layer decreases the superfluid current of e — h 
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pairs in the bilayer system. This decrease is caused by correlations between 
spatial fluctuations of superfluid density n s and superfluid velocity v s . namely, 
under influence of the donor electric field the superfluid density n s becomes in- 
homogeneous, and due to the current conservation law divj s = div(n s v s ) = 
in the regions where n s is higher the value of v s is lower. These correlations 
lead to decreasing of the average superfluid current density in the second order 
of the perturbation theory. 

The following analysis is based on the nonlinear equation for the wave func- 
tion of the e — h pair condensate obtained in Ref. 10: 

..d$(R,t) tf d a(R)r i\ 1 / d a(R) r i\ 

-ia(R)i?L + ^(R)-^ + .9(R)l*(R,i)| 2 }*(R,0- (1) 

Note that the equation (p} is obtained within the self-consistent field approx- 
imation at low density of e — ft. pairs (yx -C 1). The equation ([T]) is similar to 
the Gross - Pitaevskii equation [13], that is not surprising since at low density 
the e — h pairs act as a weakly nonideal Bose gas [UMTS]. Due to inhomogeneity 
of the bilayer system the coefficients in (pj (the effective mass of the pair Mh, 
its polarizability a and the interaction constant g) are functions of coordinates 
in the plane of the layers: R = (x,y). Explicit expressions for these coeffi- 
cients and for the potential energy of the e — h pairs U (R) will be shown below. 
The significant difference of ([T]) from the Gross - Pitaevskii equation consists 
in changing the structure of the first term in the right hand side. Firstly, the 
action of the gradient operator spreads to the coordinate dependent mass of 
the e — h pair, and secondly, the gradient operator itself is "extended" with an 
addition related to the polarizability of the e — h pairs. This addition contains 
the effective electric field E tof = E|| + Ej n t, where Em is the electric field of the 
external charges (parallel to the layers) and Ei„ t is the field caused by curvature 
of the conducting layers [TUl [16] . 

It is easy to demonstrate that the imaginary part of the equation ([lj is re- 
duced to the continuity equation n s -|-divj s = 0, where the superfluid component 
density is n s =| 'J | 2 and the superfluid current density is 

Here <p is the phase of the wave function =| ^> | e IV ). 

Now we consider the case when the inhomogeneity is caused by random 
positions of ionized donor atoms in a layer at a distance of Zd from the bilayer 
system (Fig. [1]). In such a situation the coefficients in Eq. ([1} do not depend 
on coordinates and Etot coincides with the parallel electric field E||. The field 
En and the potential energy of an e — h pair £/(R) are expressed in terms of 
the electrostatic potentials Vi(R) and V^(R) created by the charged donors in 
the electron and hole layers correspondingly: Ey = — (1/2) ^-[Vi (R) + V2(R)]; 
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Figure 1: Schematic view of a bilayer electron- hole system with interlayer dis- 
tance d. Electron-hole pairs with spatially separated components and a layer of 
donor atoms located at a distance Zd from the bilayer system are shown. 

C/(R) = e[V2(R) — Vi(R)]. In this case the real part of the equation (fTJ) can be 
written in the form 



2M H s/rTs <9R 2 2 s 2 



V US ~ + (R) 2 + e[V 2 (R) - Vt(R)] - gn s + » = 0. (3) 



Note that experiments with bilayer systems are carried out in strong mag- 
netic fields, when the magnetic length is much less than the Bohr radius of the 
e — h pair. In this limit the effective mass of the e — h pair Mh appears as the 
result of Coulomb interaction between the electron and the hole [T71 ITS] and its 
value does not depend on the intrinsic masses of the carriers: 



M. 



H 



2eh 2 

C 2 Ih V 7T 



( 1 + f) eXP (2%) erfC ( 



/ d 2 



2 di- 1 

7T l H . 



Here e is the dielectric constant and erfc(x) = 1 — -4L f Q 



(4) 

*" dt is the com- 



plementary error function. Besides Mjj, the equation ([3]) also contains the 
polarizability of the e — h pair a = c 2 Mh/H 2 and the interaction constant 



9 = 



^d-(2Kf' 2 l H + {2^f' 2 l H exp(^)erfc(-^-) 



(5) 



The expression for the electrostatic potentials created by donors in the elec- 
tron and hole layers has the form 



V m = - / dR'p(R') (R-R') 



(6) 
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where the density of donor atoms is p(R') = Yli 5(R' — R-)- In this equality R^ 
determines the position of the i-th donor. 

If the potentials Vi/ 2 )(R) change signihcantly at distances much larger than 
the coherence length £ = H/ 1 \JMngn s then the first term in (|3]) can be neglected 
as small relatively to gn s . Since we are interested in the linear in Vj.(2) correction 
to n s , we also omit the second and third terms which are quadratic in n s . As a 
result we obtain the equation for n s in the Thomas - Fermi approximation: 

e[V 2 (R) - Vi(R)] - gn s (R) + /i = 0. (7) 

It is convenient to split the potential difference V\ (R) — V% (R) to the homo- 
geneous part (Vi — V-z) = j* J c?R[Vi(R) — V2(R)], where S is the system area, 
and the deviations from the average value V\ — V2 = Vx — V%— (Vi — V2). The 
homogeneous part determines the average filling factors v\ and v-i , whereas the 
inhomogeneous part gives the spatial fluctuations of the superfluid density of 
e — h pairs we are interested in. 

n,(R) = -[Vi(R)-Va(R)]. (8) 
9 

Since n s and v s are related by the equation div(n s v s ) = 0, the density fluc- 
tuations fi s lead to fluctuations of the superfluid velocity v s . In the linear 
approximation we have 

Oti 

v «o-7jj^ +n s0 divv s = 0, (9) 

where the index "0" denotes the homogeneous part. This scalar equation is 
insufficient to find v s . An additional equation can be obtained from @. In the 
absence of quantized vortices 

Cllll * s = -^ dW ^- (10) 
In order to find v s we calculate curls of both sides of this equality. For the 
Fourier component of the superfluid velocity w s k we obtain 

fc 2 v sk - k(kv sk ) = ^_(kE|| k )[kH]. (11) 
cMh 

Equation ([9]) leads to (kv s k) = — (kv s o)'«sk/"-sO- Substituting this expression 
into (fTTj) yields 

n sk k(kv s0 ) a (kEjk) 

Thus, fluctuations of the superfluid velocity v s are determined by fluctuations 
of the superfluid density n s and fluctuations of the electric field E||. 

Since in real bilayer systems Zd 3> d, when calculating the Fourier compo- 
nent of the parallel electric field Eyk we may assume d = 0. Using the definition 
of Eyk, we obtain from ([6]): 

E,,k = — Pke-^r, (13) 
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where pk = e lkR i jf we expand Vim by d, the Fourier component of the 
superfluid density can be easily found from (j8]): 

n sk = Pke fcZd . 14 

After substituting the formulae obtained above into the expression for the 
average superfluid current density (j s ) = (n s v s ) = n s0 v s0 + Sk/o n skV s -k we 
obtain 

0.) = nsov s0 - J- E(— ) Vp- k e-^M^) + 
' k^o y 

E47r 2 e 2 da [kH] . „, 

k#0 e 2 gcM H k 

Here and below it is assumed that the area of the system 5=1. After averaging 
(|15p over positions of the donor atoms (which is reduced to replacing pk/°-k with 
the donor density rid) and changing summation by k to integration, we note that 
the second integral becomes zero due to integration over angles of the vector k. 
Calculating of the remaining integral does not present difficulties and yields 

7r / dc^ \ 2 rid 1 1 
4 V eg ) n s0 J 

The expression in square brackets is the effective superfluid density of e — h 
pairs, n e f. As follows from (|16[) . n e t < n S Q, thus, under the action of the 
electric field of randomly positioned donors the average density of superfluid 
current decreases. 

Now we consider in detail the case of small interlayer distances. At d <§; Ih 
the interaction constant of the e — h pairs 

g = W^, (17) 

el H 

and for the effective density of the superfluid component we find an expression 

Since rid/n s o ~ 1 and Ih 3> d. in this case a superfluid state exists 

only at quite large distance of donor layer from the bilayer system. Another 
aspect of the obtained result is that at given Zd the destructive influence of 
donors can be reduced by decreasing Ih to values about d and increasing v\ 
to the maximum value 1/2. Note that just at d ~ Ih and v\ = v% = 1/2 
bilayer systems possess properties which can be explained by appearance of 
superfluidity of e — h pairs El [5] ■ 

Since two layers of donor atoms are used in experiments, one of which is 
locates above the bilayer system and another is under it, the expression (|16|) 
must be generalized to this case. Calculations (analogous to that for one layer) 
yield 



G 



{is) = 



v s0 , (19) 



where n<n is the density of donor atoms in the i-th layer, Z,n is the distance 
from the bilayer system to the i-th donor layer (i = 1,2). 

3 Suppression of superfluidity of e — h pairs by 
variations of interlayer distance 

In this section we will show that in bilayer systems with curved layers variations 
of the interlayer distance lead to a decrease of the average superfluid current 
density. So, let the interlayer distance be a smooth function of coordinates in 
the plane of the layers: d = d(R). Since the effective mass of a pair M#, its 
polarizability a and the interaction constant g depend on d, it turns out that 
they also depend on R. In the absence of electric fields created by external 
charges the potential energy of the e — h pair coincides with its binding energy 
which, according to [TO], has the form 

Eo =-£/l exp (4) erfc (^)- ^ 

Owing to the dependence of Eq on d a random potential profile occurs in the 
inhomogeneous system: £/(R) = Eo[d(R)]. Finally we remind that curvature of 
the conducting layers leads to a field ~E int which polarizes the e — h pair [TU] : 

Em * = 2k^ F{d) M [zi{K) + Z2(R)] - (21) 

Here z 1 ( 2 )(R) denotes the z coordinate of the electron (hole) layer, and the 
dimensionless function F(d) is determined by the expression 



W-4[l-V5^-p(^)-*(^)]. (22) 

In the inhomogeneous bilayer system where the effective mass of a pair is 
a function of coordinates it is convenient to write down the expression for the 
superfluid current density as j s (R) = i , C(R)p s (R). Here the coefficient 

« R > - mitr <23 » 

and the kinematic momentum of the e—h pair 

P*(*) = 4g + ^M4 (24) 
Oti, c L J 

To find (j s (R)) we perform calculations analogous to the calculations in 

the previous section. Zu2) (R) can be expressed as a sum of the average value 

and the fluctuation addition: zu2)(R) = ( z i(2)) + -Zi(2)(R)- I n the linear order 

in fluctuations the current conservation condition divj s = gives a relation 

between Fourier components K k and p sk : 
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(kp sk ) = -(kp s0 )A" k /-^o- 
Calculating now curl curl p s and using (f2"5"j) , we obtain 

(kp s0 ) , 



Psk 



EintkH 



(25) 



(26) 



K 

where the polarizability a is taken in the zeroth order in fluctuations. After 
substituting p s _k into the expression for the average superfluid current (j s ) = 
Kop s o + J2k^o ^kPs -k, we have 



An fc C ' 

k^O k^O 



E, 



it -k 



H 



(27) 



To find the superfluid density n s appearing in the coefficient K (|23[) . one 
must use the equation ([1} . For the case considered the Thomas - Fermi approx- 
imation reduces the real part of ([T]) to the equality 



£ [d(R)]+.g(R)n s (R)- M = 0. 



(28) 



Equation ([28| can be split to the homogeneous and fluctuating parts. The 
homogeneous part determines ro s n, and the fluctuating part gives h s - deviations 
of the superfluid density from to the average value. In the linear approximation 
in fluctuations 



n.(R) 



UdEo 
9 



n S Q- 



dg 



dd " v dd 

where d(R) = Zi(R) — z-z (R). From (f2"3")l we obtain 



d(R), 



A k = - 



™ S 



M 



f 



3 5 



f 8Mh 



Mm dd 



The expression for Ei„ t k directly follows from (J2T 



(29) 



(30) 



E 



mi k 



2el 



H 



-F(d)ik(zik + z 2 k). 



(31) 



After substitution of ([30]) and (|3"Tj) into (j s ) integration over angles of the 
vector k reduces the last term in (|27l) to zero. The remaining terms yield 



Thus, variations of the interlayer distance lead to a decrease of the average 
current. It follows from (|32p that in a bilayer system with curved conducting 
layers the effective superfluid density is 

r lr 1 /03d 3<A , 1 dM^ 2 , 
^ = „ rt |l--[_^_ +rh0 _j + __j <d)j. (33) 

At d <C this expression is significantly simplified: 
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M V~^?~di~)- (34) 



Comparison of the expressions (JT5J) and (IM)) gives that variations of the 
interlayer distance suppress the superfiuid density stronger than ionized donors, 
if 

y > it ^ 

If we take typical values = 300 nm, Ih = 30 nm, d = 20 nm, we find (d 2 ) 1 / 2 > 
1 nm. Thus, in the case d < Ih we can assume that the influence of interlayer 
distance fluctuations prevails over the influence of ionized donors. Being based 
on this estimation, in the next section we will consider the suppression of the 
superfiuid transition temperature caused by variations of d. Note that within 
the approach developed here simultaneous consideration of influences of ionized 
donors and interlayer distance fluctuations in T c is also possible, although it 
leads to quite cumbersome expressions. 

4 Influence of interlayer distance fluctuations on 
the superfiuid transition temperature 

It is well known that in a 2D system of neutral bosons superfluidity is destroyed 
when formation of vortices leads to decrease of the free energy (Berezinskii - 
Kosterlitz - Thouless transition) [19l [20]. The transition temperature T c is 
found from the equation 

T c = ^~K(T C ). (36) 
2k B 

Since e — h pairs at low density behave as a gas of Bose particles, the equality 
(|36p also determines the superfiuid transition temperature in the bilayer systems 
considered here. It is important to note that the value of K decreases both due 
to collective excitations of the condensate (phonons) and under the action of 
inhomogeneities. Calculation of the influence of phonons does not cause diffi- 
culties (see details e. g. in [21]), and the influence of d(R) has been considered 
in the previous section. Combining both contributions, we obtain 

3C(3)fc|T 3 1 (dK^ ~ 
K{T) ~ K ° 2^Vn s0 " 2K Q \~dd) {d } ' (37) 



where <(3) = 1,202 and 



dK 



1 fdE dg\ 1 dM H 

I 1- n „r. -I 



gn s a V dd ' s ° 'ddJ Mh dd 



(38) 



Substituting K(T) into (|36j) yields a cubic equation for the critical temper- 
ature: 



9 



Here 



3 _ 2Trh 2 g 2 n 2 s0 
1 o — 



3C(3)A;| 

If the following notation is used: 



3nh 2 



K »-2^(f)>> 



2A: 

the critical temperature T c can be represented by a compact formula 



T c = 7} 




1/^3 
4 VT* 



1/2 




4 V: 



It is convenient to analyze the results in the graphic form. First we assume 
d = and find the dependence T C (H, d) for a homogeneous bilayer system (Fig. 
[3]) . It follows from Fig. 2 that at a given value of the magnetic field the critical 
temperature as a function of the interlayer distance has the form of a curve 
with a maximum. The maximum value of T c is reached at d = d ma x, and 
the dependence d max (H) can be approximated quite well by the curve d max = 
0, <Uh- In the d < d max region the critical temperature rapidly decreases due to 
strengthening thermal fluctuations with decreasing the interaction constant g. A 
smooth decrease of the critical temperature at d > d max is caused by increase of 
the effective mass of the e — h pair as d increases. Note that this result describes 
the behavior of T c (d) only qualitatively, because at large d destruction of phase 
coherence is first of all caused by quantum fluctuations |11| . i. e. due to an 
effect which is beyond the self-consistent field approximation used by us. 

Now we proceed to consideration of bilayer systems with random variations 
of the interlayer distance d ^ 0. The influence of such variations on the su- 
perfluid transition temperature is shown in Fig. 3, where the curves of T c (d) 
are shown at different values of the mean square fluctuation (d 2 ) 1 / 2 . It is obvi- 
ous that random inhomogeneities of the interlayer distance suppress the critical 
temperature. This suppression has such a character that the form of the T c (d) 
curves is unchanged, but the height of the maximum of T c (d) significantly de- 
creases and its position shifts to the range of larger values of d. While analyzing 
Fig. ??, of course, one must take into account that the behavior of T c in the 
range, where the corrections to the critical temperature are of the order of T c 
itself, must be realized only as qualitatively correct. 

If we consider each curve in Fig. 3 as a boundary between superfluid and 
normal phases at corresponding {d 2 ) 1 / 2 . we can make a conclusion that at given 
values of temperature and magnetic field the superfluid phase exists in a par- 
ticular range of interlayer distances, d c \ < d < d C 2, and the size of this interval 
decreases with increasing the mean square fluctuation [d 2 ) 1 / 2 . It is obvious 
from Fig. 3 that the decrease of the superfluidity existence range is mainly 
caused by strong influence of inhomogeneities on d c \. 
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Figure 2: Temperature of the superfluid transition T c as a function of magnetic 
field H and interlayer distance d. The temperature T c is calculated for a homo- 
geneous system at fixed filling factor v\ = 0.1 and dielectric constant s — 12.5. 
The domain of applicability of the result is located to the left and lower than the 
dashed line given by equality d = Ir\ the dash-dot line d — OAIh approximates 
the positions of maxima of the dependence T c (d) at a given magnetic field (see 
main text). 
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Figure 3: Influence of the interlayer distance on the superfluid transition tem- 
perature. Curves T c (d) are calculated for a magnetic field H = 5 T and mean 
square deviations (d 2 ) 1 / 2 = 0, 0.3, 1 and 2 nm (top-down). 
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5 Summary 



The influence of the randomly distributed inhomogeneities on the superfluidity 
of the e — h pairs in bilayer systems has been considered in this article. An 
important and partly unexpected result is a high sensitivity of systems with 
small interlayer distances (d <C Ih) to inhomogeneities. Formally, this is ex- 
pressed through the appearance of the factors (d/ln) 2 in expressions (fT5| and 
(l3~4")l which describe the decrease of the effective superfluid density n e f due to 
the electrostatic field of the ionized donors and the roughening of the conducting 
layers respectively. The physical origin of the strong suppression of n e f is the 
weakness of the interaction between the e — h pairs for d <C Z# (see Eq. pT]) ). 

The results of this article have been obtained under assumption that the 
filling factor of the ground Landau level v\ <C 1. From our point of view 
this limitation is not significant, since according to equalities (|40|) and (|4"Tj) the 
increase of the pair density n s o (i. e. the increase of v\) leads to the increase of 
the critical temperature T c , but does not change the picture represented in Fig. 
3. Thus, one can expect that these results qualitatively correctly describe the 
behavior of T c (d) in the inhomogeneous systems with filling factors v\ k, 1/2 
which are commonly used in experiments on finding the superfluidity of the e — h 
pairs. 

Also we note that several important conclusions for future experiments can 
be drawn. Firstly, one of the critical conditions for observing superfluidity in 
bilayer electron-hole systems is a high degree of homogeneity of these systems. 
Secondly, in homogeneous systems the maximum value of T c is reached for the 
interlayer distance d « OAIh- The superfluid phase of the interwell excitons has 
to be looked for in the vicinity of this value. Additionally one has to keep in 
mind that in the real, i. e. inhomogeneous, systems the range of the interlayer 
distances d at which the observation of the superfluid phase is most probable, 
is shifted towards higher values (see Fig. 3). 
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